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1)  ABSfnACf  (AMiuMi/fM  ?Q0  ivof  di) 

A  "Workshop  on  Problems  of  Rotating  Fluids"  was  held  at  the  U.S.  Army  High 
Performance  Computing  Research  Center  on  22-23  of  April  1991.  This  meeting 
was  co-sponsored  by  M.  Miller,  CREEC  and  D.  Joseph,  University  of  Minnesota. 
The  purpose  of  the  workshop  was  to  review  the  status  of  technology  related 
to  rotating  fluid  dynamics  with  special  emphasis  on  liquid-filled  projectiles. 
Sixteen  technical  papers  were  presented  by  individuals  frcm  the  government  and 
academia . 
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FORWARD 


The  "Workshop  on  Problems  of  Rotating  Liquids"  was  held  at  the 
Army  High  Performance  Computing  Research  Center  (AHPCRC)  in 
Minneapolis,  MN  on  22-23  April  1991.  The  two  day  meeting  was 
co-sponsored  by  Professor  Daniel  Joseph,  University  of  Minnesota 
and  Miles  Miller,  U.S.  Army  Chemical  Research,  Development  and 
Engineering  Center  (CRDEC) .  Dr.  Joseph  is  engaged  in  fluid 
dynamic  research  under  an  Army  Research  Office  (ARO)  contract 
and  Mr.  Miller  is  the  Scientific  Coordinator  tor  the  Fluid 
Dynamics  Area  of  the  CRDEC  Basic  Research  Program. 

The  "Workshop"  primarily  dealt  with  the  subject  of  liquid-filled 
projectile  flight  1.. suabilities.  The  meeting  location  was 
selected  to  reflect  the  critical  use  of  Computational  Fluid 
Dynamics  (CFD'  and  high  performance  computers  in  support  of 
this  technology  area  and  was  the  first  "Workshop"  hosted  by 
the  AHPCRC.  Those  invited  to  attend  had  all  contributed  to 
the  technical  aspects  of  this  topic.  The  list  of  attendees 
is  contained  on  page  337. 

The  two  general  goals  of  the  "Workshop"  were  to: 

*  Help  resolve  an  important  problem  in  the 
nation's  interest. 

*  Move  the  general  subject  of  fluid  dynamics 
forward . 

More  specific  objectives  of  the  "Workshop"  were  to  assess  the 
progress  made  in  this  technology  and  recommend  research 
activities  to  address  particular  needs. 

The  "Workshop"  agenda  is  shown  on  pages  339  and  340.  The  first 
day  was  devoted  to  overviews,  tours,  and  special  presentations. 
The  second  day's  activities  were  concerned  with  various  "work 
in  progress"  sponsored  or  performed  by  either  the  U.S.  Army 
Ballistic  Research  Laboratory  (BRL)  or  the  CRDEC.  The  papers 
presented  by  individuals  from  the  Ohio  State  University 
represent  work  performed  under  a  research  contract  funded  by 
the  CRDEC.  These  proceedings  contain  copies  of  the  viewgraphs 
shown  by  each  presenter. 

The  last  time  a  special  meeting  on  this  subject  was  held  was 
at  the  "Roundtable  On  Liquid-Filled  Shell"  held  in  September  of 
1984.  That  meeting  had  charted  the  course  of  study  in  this 
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area  for  the  past  seven  years  and  It  seemed  timely  that  we 
should  reassess  our  progress  at  this  point.  During  the 
"Wrap-Up"  session  of  the  "Workshop",  the  recommendations  of  the 
previous  "Roundtable"  were  reviewed  by  the  attendees  and  an 
updated  list  of  new  topics  and  future  directions  for  research 
in  this  area  was  formulated.  These  results  are  summarized  in 
the  last  section  of  these  proceedings. 
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PREFACE 


The  Workshop  on  Problems  of  Rotating  Liquids  was  held 
22-23  April  19S1  at  the  U.S.  Army  High  Performance  Computing 
Research  Center  (AHPCRC),  Minneapolis,  MN,  and  was  authorized 
under  project  No.  FI -1-1Z 1RA-NMCL. 

The  use  of  trade  names  or  manufacturers'  names  in  this 
report  does  not  constitute  an  official  endorsement  of  any  commercial 

*  products.  This  report  may  not  be  cited  for  purposes  of 
advert i sement . 

Reproduction  of  this  document  in  whole  or  in  part  is 

•  prohibited  except  with  permission  of  the  Commander,  T'.S.  Army 
Chemical  Research,  Developmert  and  Engineering  Center, 

ATTN:  SMCCR-SPS-T ,  Aberdeen  Proving  Ground,  MD  21010-5423. 

However,  the  Defense  Technical  Information  Center  and  the  National 
Technical  Information  Service  are  authorized  to  reproduce  the 
document  for  U.S.  Government  purposes. 

Thi;  report  has  been  approved  for  release  to  the  public. 


Acknowledgments 

Special  thanks  is  given  to  Dr.  Sell  and  his  staff  at 
the  AHPCRC  for  making  the  various  lodging  and  meeting  arrangements. 


1  Accession  For  j 

NT  IS 

or  aai 

DT1C 

T.\!3 

n 

LV'ii.-r 

cu:ic--d 

□ 

Ji.it  1 

f 1 4 

on _ _ 

Hv 

Dim  r 

O  U* 

Ay,-,  l 

i.  <tL  I  j.  i 

■>  Cov’f.S 

* ; 

.  .  ■ 

Dlst 

;  :\l 

_ 

t 


Bl  ank 


6 


CONTENTS 


Page 


OVERVIEW  OF  CRDEC  RESEARCH  PROGRAM  FOR  LIQUID-FILLED 

PROJECTILES  .  9 

ANALYSIS  AND  VISUALIZATION  OF  THE  FLOW  IN  A  SPINNING 

AND  NUTATING  CYLINDER  .  53 

A  CENTRIFUGAL  SPECTROMETER  .  121 

MOMENT  EXERTED  BY  A  VISCOUS  LIQUID  IN  A  SPINNING, 

CONING  CONTAINER  .  131 

NUMERICAL  STUDY  OF  UNSTEADY  3D  FLOWS  IN  A  SPINNING 

AND  NUTATTNG  CYLINDER  .  159 

NUMERICAL  SIMULATION  FOR  NON-CYLI NDRICAL  LIQUID-FILLED 

CONTAINERS  .  175 

MOTION  OF  TWO  IMMISCIBLE  FLUIDS  IN  A  SPINNING  AND 

CONING  CYLINDER  .  199 

DIRECT  MEASUREMENT  OF  LIQUID  EFFECTS  USING  A  MOMENT 

BALANCE  .  231 

LABORATORY  FLIGHT  STABILITY  EVALUATION  OF  PRODUCTION 

M825A1  PAYLOAD  CANISTERS  .  251 

THEORY  AND  EXPERIMENTS  FOR  ROTATING  POROUS  MEDIA  FLOW  ....  271 

EFFECTS  OF  INTERIOR  CANISTER  WALL  ROUGHNESS  ON  LIQUID 

DESPIN  MOMENTS  .  295 

INSTRUMENTED  FLIGHT  TESTS  ARTILLERY  PROJECTILES  WITH 

SELECTED  LIQUID-FILLS  .  309 

SUMMARY  OF  NEW  DIRECTIONS  FOR  LIQUID-FILLED  PROJECTILE 

STUDIES  .  333 

"WORKSHOP"  RECOMMENDATIONS  (Not  Prioritized)  .  335 

LIST  OF  REGISTRANTS  .  337 

PROGRAM  AGENDA  .  339 


7 


* 


Bl  ank 


8 


ROTATING  FLUIDS  WORKSHOP 


X  O 
O  uj 
cc 

<  o 

LU  CC 

(/>  Q. 

g  Q 

u-  |jj 

O  “i 

LU  =f 
Q  Li. 


Sg 

>  Li. 


5  <£ 

?  « 
.S  a> 

O)  g 
c 

LU  £ 


«J  £  2 

£  ®  m 
a>  E  0) 

£  V  “ 

“■2  m 
0) 

>**■</> 
o>  ^  <5 

*.  Q  w 

■g  S5  a 
k-  .r  a> 
«o  co  o 
a>  >*  c 
&  jz  o 

«  a  o 

«_  w  -D 

8  5  g 
E  *g 

<P  V- 
JC  <0 
O 

w  w 

£c 


Army  High  Performance  Computing  Research  Center 

University  of  Minnesota 


OVERVIEW  OF  CRDEC  RESEARCH  PROGRAM 
FOR  LIQUID-FILLED  PROJECTILES 


In  the  short  time  I  have  this  morning,  I  want  to  indicate 
the  reasons  why  the  CRDEC  was  interested  in  holding  this 
workshop;  what  oar  main  technical  interests  are;  and  provide 
a  brief  review  of  the  background,  current  efforts,  and 
future  directions  in  this  area. 

The  main  area  of  "Rotating  Liquids"  of  interest  to  the  CRDEC 
has  to  do  with  the  flight  stability  of  liquid-filled 
projectiles.  The  CRDEC  is  the  development  agency  for 
chemical  weapons  which  include  antipersonnel,  antimateriel, 

flame  &  incendiary  and  smoke  &  obscurants.  Most  of  these  » 

weapons  involve  spinning  projectiles  which  contain  liquid  or 
non-rigid  payloads.  This  would  include  artillery  shells, 
mortars,  grenades,  small  arms,  missiles,  etc. 

At  the  previous  "Roundtable  on  Liquid-Filled  Projectiles", 
held  in  September  1984,  a  relatively  detailed  history  of 
this  technology  was  presented.  For  the  benefit  of  those  of 
you  who  are  new  to  this  topic,  I  would  like  to  briefly 
review  certain  highlights. 

This  viewgraph  depicts  a  historical  perspective  of  the 
significant  events  in  this  area.  While  the  CRDEC  has  been 
interested  in  this  problem  since  World  War  I,  until  the  late 
1970 | s ,  their  effort  was  limited  primarily  to  experimental 
studies  in  direct  support  of  developmental  munitions.  The 
bulk  of  the  scientific  research  was  accomplished  by  the  BRL. 

The  Stewartson-Wedemeyer  Theory  appeared  to  adequately 
address  the  liquid-fills  of  interest  at  that  time. 

Consequently,  the  BRL  effort  concentrated  on  experimental 
activities  such  as  free  gyroscope  and  yaw  sonde 
investigations . 

The  research  activities  at  the  CRDEC  started  in  1977  as  a 
consequence  of  a  flight  stability  problem  with  XM761;  a 
155mm,  artillery  projectile  which  had  a  partial 
solid/partial  liquid  payload  composition.  The  result  of 
this  effort  was  the  evolution  of  the  Test  Fixture  For 

Non-Rigid  Payloads.  This  experimental  apparatus  forces  a  m 

full  sized  payload  container  to  assume  the  simultaneous 

spinning  and  coning  motion  of  a  projectile  in  flight.  The 

test  fixture  was  used  to  eliminate  the  XM761  problem 

culminating  in  the  M825  projectile.  It  was  subsequently 

employed  to  discover  and  characterize  the  destabilizing  * 

influence  of  highly  viscous  liquid-fills. 
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The  real  beginnings  of  the  CRDEC  research  effort  took  place 
in  1982  with  the  establishment  of  the  Fluid  Dynamics  Work 
Area  of  the  CRDEC  Basic  Research  Program.  At  the  same  time, 
the  Laboratory  Test  Fixture  For  Non-Rigid  Payloads  was 
extensively  modified  to  increase  its  performance,  automate 
the  data  acquisition  and  reduction,  and  facilitate 
experimental  operations.  It  is  currently  being  employed  for 
basic  research  studies  as  well  as  for  various  developmental 
munition  programs.  Examples  of  both  uses  will  be  presented 
as  part  of  the  "Work  In  Progress"  portion  of  this  meeting. 

The  significant  point  of  this  viewgraph  is  the  extensive 
amount  of  work  put  into  this  technical  area  by  the  army 
during  the  last  ten  years  and  the  complimentary  efforts  of 
the  many  different  individuals  and  organizations  involved. 
Their  combined  efforts  have  caused  the  technology  associated 
with  this  problem  to  progress  considerably  as  illustrated  in 
the  next  viewgraphs. 

The  destabilizing  moment  produced  by  a  Liquid-fill  in  a 
spinning  and  coning  cylinder  can  be  depicted  as  a  function 
of  the  fluid  characteristics  (i.e.,  Reynolds  number)  and  the 
projectile  motion  (i.e.,  non-dimensional  coning  rate). 

Fifteen  years  ago,  the  theory,  and  in  fact  the  knowledge,  of 
liquid-filled  projectiles  was  limited  to  very  low  viscosity 
liquid-fills.  These  could  be  predicted  quite  nicely  by  the 
Stewartson-Wedemeyer  Theory.  Ten  years  ago,  we  were  aware, 
through  laboratory  experiments,  cf  a  new  form  of  instability 
at  the  higher  viscosities,  but  had  no  predictive  theory. 
Within  the  last  decade,  we  have  developed  theoretical 
methods  to  handle  the  entire  range  of  liquid  viscosities. 
Murphy  has  extended  the  Stewartson-Wedemeyer  theory  for  the 
high  Reynolds  numbers  and  computerized  its  use.  Hall, 

Sedney  and  Gerber  and  Herbert  and  Li  have  both  developed 
theoretical  methods  which  cover  the  lew  Reynolds  number 
range.  The  general  capabilities  of  the  various  analytical 
methods  are  summarized  in  the  next  viewgraph. 

The  CRDEC  Basic  Research  Program  has  followed  a  three 
pronged  approach  as  illustrated  in  this  viewgraph.  The 
theoretical  and  numerical  studies  are  completed  under 
research  contracts  through  academic  institutions  with  the 
CRDEC  performing  primarily  experimental  investigations  using 
their  Laboratory  Test  Fixture  for  Non-Rigid  Payloads.  The 
following  series  of  viewgraphs  summarizes  the  significant 
accomplishments  achieved  through  the  CRDEC  sponsored 
program.  The  next  viewgraphs  list  the  publications  related 
to  this  area  which  were  sponsored  by  the  CRDEC.  This  work 
has  been  complimented  by  research  at  the  BRL  and  various  ARO 
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funded  research.  CRDEC  and  BRL  personnel  meet  every  two 
months  to  review  their  respective  activities. 

The  CRDEC  effort  has  concentrated  on,  but  not  been  limited 
to,  highly  viscous  liquids.  The  goal  has  been  to 
understand,  predict  and  prevent  flight  instabilities  due  to 
any  liquid-fill.  A  comprehensive  summary  of  papers  on  the 
year's  research  findings  are  presented  each  year  at  the 
CRDEC  Scientific  Conference  on  Chemical  Defense  Research.  A 
review  of  the  Fluid  dynamics  progress  is  held  for  the  AK'' 
each  spring  at  the  CRDEC.  The  following  is  an  abbreviated 
version  of  this  year's  CRDEC/ARO  review  which  took  place  in 
early  April  1991. 

The  overall  objective  of  the  Fluid  Dynamics  Work  Area  of  the 
CRDEC  Basic  Research  Program  is  shown  in  this  viewgraph. 
While  the  area  of  antipersonnel,  chemical  munitions  has  been 
reduced,  it  has  been  replaced  with  the  other  areas  as  shown. 
While  efforts  still  continue  to  understand  and  predict 
flight  instabilities,  the  current  emphasis  is  on  determining 
ways  to  reduce  or  eliminate  them. 

The  Fluid  Dynamics  Area  involves  two  areas  of  emphasis: 
Fluid-Filled  Projectile  Flight  Instabilities  and  Fluid 
Rheology.  The  latter  topic  is  important  not  only  for  this 
area,  but  for  other  CRDEC  programs  as  well .  The 
Fluid-Filled  Projectile  work  receives  the  major  portion  of 
the  funding  for  the  reasons  shown  in  this  viewgraph:  problem 
unique  to  CRDEC;  CRDEC  possesses  special  experimental 
facilities;  and  new  chemical  munition  fills  must  be 
addressed.  As  indicated  in  this  viewgraph,  the  CRDEC 
in-house  effort  deals  mainly  with  experimental  studies  which 
support  theoretical  studies  performed  under  contract  by 
universities 

This  viewgraph  contains  selected  accomplishments  during  FY91 
in  the  theoretical/numerical  area  achieved  by  Dr.  Herbert, 
Dr.  Li  and  Mr.  Selmi,  researchers  at  the  Ohio  State 
University  under  a  CRDEC  research  contract.  The  next  two 
viewgraphs  highlight  the  computer  graphics  techniques  and 
the  two,  immiscible  fluid  analysis  being  pursued.  In 
particular,  the  two  fluid  analysis  includes  the  case  of  very 
small  amounts  of  the  lower  viscosity  fluid  which  is  the 
basis  for  the  additive  approach  to  reducing  the 
destabilizing  moment  due  to  the  higher  viscosity  fluid. 

FY91  experimental  accomplishments  at  the  CRDEC  are 
summarized  in  the  next  viewgraph.  The  next  three  viewgraphs 
include  additional  details  for  these  items.  Detailed 
studies  were  conducted  on  the  CRDEC  test  fi  :ture  to  assess 


the  effect  of  relative  viscosity,  density  and  amount  of 
additive.  These  data  can  be  compared  with  the  theoretical 
results  shown  previously.  Test  Fixture  experiments  showed 
that  the  payload  container  inner  surface  roughness  has  no 
influence  on  the  flight  stability.  The  CRDEC  Test  Fixture 
was  used  to  evaluate  an  artillery  projectile  containing  a 
viscoelastic  fill  for  a  special  Operation  Desert  Storm 
application.  The  ability  to  measure  the  despin  moment  due 
to  an  arbitrary  non-rigid  payload  on  the  Test  Fixture  and 
then  to  compute  the  associated  destabilizing  yawing  moment, 
with  confidence,  was  established  through  the  basic  research 
program  on  liquid-filled  projectiles. 

Additional  accomplishments  are  listed  in  the  next  viewgraph 
followed  by  four  viewgraphs  illustrating  the  basic  results 
of  the  instrumented  flight  test  program.  These  illustrate 
that  the  partial-fill  condition  produces  the  same 
instability  as  the  fully-filled  condition;  the  shear 
thinning,  viscoelastic  fill  provides  a  stable  flight;  and 
the  immiscible,  low  viscosity  additive  does  eliminate  flight 
instabilities,  but  transient  effects  must  be  considered. 

The  classic  "Epicyclic  Theory"  of  projectile  motion  has  been 
modified  to  include  the  effects  of  a  liquid-fill;  indicating 
the  combination  of  the  internal  payload  and  external 
aerodynamic  effects  on  the  projectile  flight  flight 
stability. 

Future  Directions  in  this  CRDEC  research  area  are  shown  in 
the  next  viewgraph.  Work  w.iJl  continue  to  assess  the  use  of 
immiscible,  low  viscosity  additives  to  reduce  or  eliminate 
viscous,  liquid-fill  flight  instabilities,  especially 
considering  launch  transient  effects.  The  potential  use  of 
longitudinal  baffles  as  a  means  of  eliminating  flight 
instabilities  will  also  be  evaluated.  Previous  Test  Fixture 
data,  depicted  in  the  next  viewgraph,  indicate  that  they  may 
be  effective  under  certain  conditions.  Two,  three  and  four 
section  baffles  will  be  investigated  on  the  Laboratory  Test 
Fixture  for  Non-Rigid  Payloads.  The  influence  of 
viscoelastic  effects  on  creating  very  large  despins  at 
relatively  low  spin  rates  as  shown  in  this  viewgraph  will 
also  be  studied.  Other  non-Newtonian  fluids  will  also  be 
investigated . 


SIGNIFICANT  EVENTS  IN  LIQUID-FILLED  PROJECTILE 

FLIGHT  STABILITY  TECHNOLOGY 
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SIGNIFICANT  ACCOMPLISHMENTS  OF  CRDEC 
RESEARCH  IN  LIQUID-FILLED  PROJECTILES 


1977 

1978 

1979 

1981 

1982 

1983 

1983 

1984 

1984 

1985 

1985 

1985 

1985 


CRDEC  Laboratory  Test  Fixture  built  and  used  to  solve  liquid- 
fill  flight  stability  problem  of  developmental  smoke  projectile. 

Discovery  of  new  type  of  flight  instabilty  created  by  highly 
viscous  liquid-fills  from  laboratory  fixture  experiments  at 
the  CRDEC. 

Extensive  experimental  investigation  of  viscous  liquid-fills 
using  CRDEC  test  fixture. 

Flow  visualization  studies  conducted  for  internal  flow  inside 
spirming/coning  cylinder  using  CRDEC  test  fixture. 

Fluid  Dynamics  Work  Area  established  in  CRDEC  Basic  Research 
Program. 

Direct  laboratory  measurement  of  the  destabilizing  yawing 
moment  due  to  a  highly  viscous  liquid-fill  using  the  CRDEC 
test  fixture. 

Computational  Fluid  Dynamic  (CFD)  analysis  of  the  viscous 
liquid-fill  situation  by  Sandia  Laboratories  supported  by 
the  CRDEC. 

"Roundtable  on  Liquid-Filled  Shell"  co-sponsored  by  CRDEC. 

Theoretical  analysis  of  viscous  liquid-fill  problem  by  Herbert 
(VPI)  including  use  of  volume  integral  to  compute  liquid- 
induced  moments  supported  by  the  CRDEC. 

Six  Degree-of-Freedom  flight  motion/trajectory  program  for 
liquid-filled  projectiles  evolved  by  Sandia  Laboratories 
supported  by  the  CRDEC. 

Finite  element  CFD  analysis  completed  for  viscous  liquid-fill 
problem  by  Rosenblat  (Fluid  Dynamics  International)  supported  by 
the  CRDEC. 

Development  of  general  relationship  between  liquid  rolling  and 
liquid  yawing  moments  by  Rosenblat  supported  by  the  CRDEC. 

Formualtion  of  a  three-dimensional  graph  to  depict  the 
entire  liquid-filled  projectile  flight  stability  problem 
by  CRDEC  researchers. 


1986  -  CFD  analysis  conducted  for  viscoealstic  fluid-filled 
projectile  by  Rosenblat  supported  by  the  CRDEC. 

1986  -  Laboratory  measurement  of  the  influence  of  a  viscoelastic 

fluid-fill  in  creating  projectile  flight  instabilities  using 
the  CRDEC  test  fixture. 

1986  -  Formulation  of  the  theoretical  relationship  between  the 

liquid-fill  induced  rolling  and  yawing  moments  by  Rosenblat 
supported  by  the  CRDEC. 

1987  -  Detailed  experiments  to  validate  the  linear  and  generalized 

liquid  moment  coeffients  using  the  CRDEC  test  fixture. 

1987  -  Fully  spectral  CFD  code  developed  by  Herbert  (Ohio  State 

Ur.iv)  for  viscous  liquid-fill  problem  supported  by  the  CRDEC. 

1988  -  Purely  analytical  method  developed  to  compute  liquid-fill 

moments  at  any  Reynolds  number  by  Herbert  and  Li  supported 
by  the  CRDEC . 

1989  -  Laboratory  measurements  of  the  effect  of  partial-fill  case 

with  a  viscous  liquid-fill  using  the  CRDEC  test  fixture. 

1989  -  Laboratory  demonstration  of  the  use  of  an  immiscible,  low 

viscosity  additive  to  eliminate  viscous  liquid-fill  flight 
instabilities.  Concept  evolved  in  collaboration  with  Joseph 
(University  of  Minnesota) . 

1990  -  Experimental  evaluation  of  effect  of  internal  surface  texture 

of  payload  container  on  creation  of  viscous  liquid-fill, 
flight  instabilities. 

1990  -  Detailed  experimental  evaluation  performed  on  influence  of 

various  immiscible,  low  viscosity  additives  using  CRDEC  test 
fixture. 

1991  -  Theoretical  analysis  completed  to  predict  effect  of  two 

immiscible  fluids  covering  entire  range  of  viscosity,  density 
and  relative  volumes  by  Selmi  (Ohio  State  Univ)  supported 
by  CRDEC. 

1991  -  Instrumented  flight  tests  conducted  to  validate  viscoelastic 
fluid,  partial-fill  case,  and  immiscible  additives. 

1991  -  Enhanced  computer  graphic  techniques  evolved  to  facilitate 
visualization  and  intrpretat  >  of  internal  fluid  dynamics 
of  spinning/coning  cyinder  b j  Herbert  supported  by  CRDEC. 


22 


BIBLIOGRAPHY  OF  LIQUID-FILLED  PROJECTILE  PUBLICATIONS 
SPONSORED  BY  CRDEC  BASIC  RESEARCH  PROGRAM 


Miller,  M.  C.;  "Flight  Instability  Test  Fixture  for  Non-Rigid 
Payloads",  Proceedings  of  1977  CSL  Technical  Conference,  Me./ 
1978. 

D'Amico,  W.  P.  and  Miller,  M.  C.;  "Flight  Instability  Produced 
by  a  Rapidly  Spinning,  Highly  Viscous  Liquid",  Journal  of 
Spacecraft  and  Rockets,  Vol  16,  No.  1,  January-February  1979. 

Miller,  M.  C.;  "Flight  Instabilities  of  Spinning  Projectiles 
Having  Non-Rigid  Payloads",  AIAA  19th  Aerospace  Sciences 
Meeting,  Paper  No.  AIAA-81-0225,  January  1981. 

Miller,  M.  C.;  "Void  Characteristics  of  Liquid-Filled  Cylinder 
Undergoing  Spinning  and  Coning  Motion",  Journal  of  Spacecraft 
and  Rockets,  Vol  18,  No.  3,  May-June  1981. 

Miller,  M.  C.;  "Flight  Instabilities  of  Spinning  Projectiles 
Having  Non-Rigid  Payloads",  Journal  of  Guidance,  Control  and 
Dynamics,  Vol  5,  No.  2,  March-April  1982. 

Miller,  M.  C.;  "Flight  Instabilities  of  Spinning  Projectiles 
Having  Non-Rigid  Payloads",  Proceedings  of  1982  Army  Science 
Conference,  June  1982. 

Miller,  M.  C. ;  "A  New  Flight  Instability  Affecting  Spinning 
Projectiles  Having  Non-Rigid  Payloads",  28th  Conference  of  Army 
Mathematicians,  Army  Research  Office  Report  No.  83-1,  February 
1983. 


Herbert,  T. ;  "Fluid  Motion  in  a  Rotating  and  Nutating 
Container",  CRDC-CR-84087 ,  July  1984. 

Miller,  M.  C.  ;  "Laboratory  Measurements  of  Despin  and  Yawing 
Moments  Produced  by  a  Viscous  Liquid  in  a  Spinning/Nutating 
Container",  Proceedings  of  the  1983  CRDC  Scientific  Conference 
on  Chemical  Defense  Research,  CRDC-SP-84014 ,  October  1984. 

Herbert,  T. ;  "The  Flow  of  a  Highly  Viscous  Fluid  in  a  Spinning 
and  Nutating  Cylinder",  Proceedings  of  the  1983  CRDC  Scientific 
Conference  on  Chemical  Defense  Research,  CRDC-SP-84014,  October 
1984. 

Vaughn,  H.  L. ,  Oberkampf,  W.  L.  and  Wolfe,  W.  P. ;  "Numerical 
Solution  to  the  Navier-Stokes  Equations  for  Spinning  and 
Nutating  Fluid  Payloads",  Proceedings  of  the  1983  CRDC 
Scientific  Conference  on  Chemical  Defense  Research, 
CRDC-SP-84014,  October  1984. 


23 


Miller,  M.  C . ;  '‘Measurement  of  Despin  and  Yawing  Moments 
Produced  by  a  Viscous  Liquid",  Journal  of  Guidance,  Control  and 
Dynamics,  Vol  6,  Ho.  2,  March-April  1985. 

Herbert,  T. ;  "Instability  of  the  Viscous  Flow  in  a  Spinning  and 
Nutating  Cylinder",  Proceedings  of  the  1984  CRDC  Scientific 
Conference  on  Chemical  Defense  Research,  CRDC-SP-85006,  June 
1985. 

Vaughn,  H.  L. ,  Wolfe,  W.  P.  and  Oberkampf,  W.  L. ;  "Six  Degree- 
of-Freedom  Simulation  of  M483  Fluid  Payload  Projectile  Using 
Numerically  Computed  Fluid  Moments",  Proceedings  of  the  1984 
CRDC  Scientific  Conference  on  Chemical  Defense  Research, 
CRDC-SP-85006,  June  1985. 

Miller,  M.  C.;  "Visualization  Studies  of  Viscous  Liquid  Flow  in 
a  Spinning  and  Coning  Cylinder",  Proceedings  of  the  1984  CRDC 
Scientific  Conference  on  Chemical  Defense  Research, 
CRDC-SP-85006,  June  1985. 

Vaughn,  H.  R.,  Oberkampf,  W.  L.  and  Wolfe,  W.  P. ;  "Fluid 
Motion  Inside  a  Spinning,  Nutating  Cylinder",  Journal  of  Fluid 
Mechanics,  Vol  150,  January  1985. 

Vaughn,  H.  R.,  Wolfe,  W.  P. ,  and  Oberkampf,  W.  L. ?  "Six 
Degree-of-Freedom  Simulation  of  Fluid  Payload  Projectiles 
Using  Numerically  Computed  Fluid  Moments",  Sandia  National 
Laboratories  Report  SAND  85-1166,  July  1985. 

Herbert,  T. ;  "Viscous  Fluid  Motion  in  a  Spinning  and  Nutating 
Cylinder",  CRDC-CR-86004 ,  November  1985. 

Miller,  M.  C.;  "Liqu  d-Filled  Projectiles  -  New  Problems,  New 
Solutions",  Proceedings  of  the  9th  American  Defense 
Preparedness  Association  International  Ballistics  Conference, 
Vol  2,  May  1986. 

Herbert,  T. ;  "Fluid  Motion  in  Liquid-Filled  Shells",  4th  Army 
Applied  Mathematics  and  Computating  Conference,  May  1986. 

Miller,  M.  c.;  "Liquid-Filled  Projectile  -  New  Problems,  New 
Solutions",  Proceedings  of  the  9th  American  Defense 
Preparedness  Association  International  Ballistics  Conference, 
Vol  2,  May  1986. 

Herbert,  T. ;  "Fluid  Motion  in  Liquid-Filled  Shells",  4th  Army 
Applied  Mathematics  and  Computing  Conference,  May  1986. 

Herbert,  T.;  "On  the  Fluid  Motion  in  Liquid-Filled  Shells", 
Proceedings  of  the  1985  CRDC  Scientific  Conference  on  Chemical 
Defense  Research,  CRDC-SP-86007 ,  June  1986. 


24 


Rosenblat,  S.;  "Viscoelastic  Fluid  Flow  in  a  Spinning  and 
Nutating  Cylinder",  Proceedings  of  the  1985  CRDC  Scientific 
Conference  on  Chemical  Defense  Research,  CRDC-SP-86007 ,  June 
1986. 

Miller,  M.  C.  and  Weber,  D.  J. ;  "Liquid-Filled  Projectile 
Flight  Instabilities,  a  Unified  Picture",  Proceedings  of  the 
1985  CRDC  Scientific  Conference  on  Chemical  Defense  Research, 
CRDC-SP-86007,  June  1986. 

Rosenblat,  S.;  "Finite  Element  Calculations  of  Viscoelastic 
Fluid  Flow  in  a  Spinning  and  Nutating  Cylinder", 

CRDEC-CR-87 02 1 ,  December  1986. 

Miller,  M.  C. ;  "Liquid-Filled  Projectiles  -  New  Problems,  New 
Solutions",  Army  R  D  and  A  Magazine,  November-December  1986. 

Herbert,  T. ;  "Viscous  Fluid  Motion  in  a  Spinning  and  Nutating 
Cylinder",  Journal  of  Fluid  Mechanics,  Vol  167,  1986. 

Miller,  M.  C.  and  Weber,  D.  J. ;  "Measurement  of  Despin  Moment 
Produced  by  a  Viscoelastic  Fluid  in  a  Spinning  and  Coning 
Cylinder",  Proceedings  of  the  1986  CRDEC  Scientific  Conference 
on  Chemical  Defense  Research,  CRDEC-SP-87008 ,  June  1987. 

Rosenblat,  S.;  "Non-Newtonian  Fluid  in  a  Spinning  and  Coning 
Cylinder",  Proceedings  of  the  1986  CRDEC  Scientific  Conference 
on  Chemical  Defense  Research,  CRDEC-SP-87008,  June  1987. 

Herbert,  T. ;  "A  Spectral  Navier-Stokes  Solver  for  the  Flow  in 
a  Spinning  and  Coning  Cylinder",  Proceedings  of  the  1986  CRDEC 
Scientific  Conference  on  Chemical  Defense  Research, 
CRDEC-SP-87008,  June  1987. 

Herbert,  T. ;  "Numerical  Study  of  the  Flow  in  a  Spinning  and 
Nutating  Cylinder",  AIAA  19th  Fluid  Dynamics,  Plasma  Dynamics 
and  Laser  Conference,  Paper  AIAA-87-1445 ,  June  1987. 

Weber,  D.  J.  and  Miller,  M.  C. ;  "Despin  Moment  Measurements 
for  Viscous  Liquids  in  a  Spinning  and  Coning  Cylinder", 
Proceedings  of  the  1987  CRDEC  Scientific  Conference  on 
Chemical  Defense  Research,  CRDEC-SP-88013 ,  June  1988. 

Herbert,  T. ;  "Flight  Simulation  for  Liquid  Filled 
Projectiles",  Proceedings  of  the  1987  CRDEC  Scientific 
Conference  on  Chemical  Defense  Research,  CRDEC-SP-88013 ,  June 
1988  . 

Miller,  M.  C.  and  Weber,  D.  J.;  "Destabilizing  Moments 
Produced  by  a  Viscoelastic  Fluid  in  a  Spinning  and  Coning 


Cylinder”,  Proceedings  of  the  1988  Army  Science  Conference, 

October  1988. 

Weber,  D.  J.  and  Miller,  M.  C.;  “Liquid  Induced  Rolling  and 
Yawing  Moment  Coefficients  for  Viscous  Newtonian  Fluids  in  a 
Spinning  and  Coning  Cylinder",  Proceedings  of  the  1988  CRDEC 
Scientific  Conference  on  Chemical  Defense  Research, 

CRDEC-SP-013,  August  1989. 

Miller,  M.  C. ;  "Laboratory  Test  Fixture  for  Non-Rigid 
Payloads",  Proceedings  of  the  International  Congress  of 
Instrumentation  in  Aerospace  Simulation  Facilities 
(ICIASF'89),  September  1989. 

Herbert,  T. ?  "Interior  Fluid  Dynamics  of  Liquid-Filled 
Projectiles",  CRDEC-CR-058,  December  1989. 

Herbert,  T.  and  Selmi,  M. ;  "Two  Fluid  Flows  in  a  Spinning  and 
Nutating  Cylinder",  Proceedings  of  the  1989  CRDEC  Scientific 
Conference  on  Chemical  Defense  Research,  CRDEC-SP-024 ,  August 
1990. 

Weber,  D.  J.  and  Miller,  M.  C. ;  "Effect  of  Partial-Fill 
Condition  on  the  Destabilizing  Moment  Produced  by  a  Viscous 
Liquid  in  a  Spinning  and  Coning  Cylinder",  Proceedings  of  the 

1989  CRDEC  scientific  Conference  on  Chemical  Defense  Research, 
CRDEC-SP-024,  August  1990. 

Miller,  M.  C. ;  "Reduction  of  Viscous  Liquid-Fill  Flight 
Instability  by  Addition  of  Immiscible,  Low  Viscosity  Liquid", 

Proceedings  of  the  1989  CRDEC  Scientific  Conference  on 
Chemical  Defense  Research,  CRDEC-SP-024,  August  1990. 

Li,  R.  and  Herbert,  T. ;  "Computational  Study  of  the  Flow  in 
Spinning  and  Nutating  Cylinders",  AIAA  Journal,  Vol  28,  No.  9, 

September  199C. 

Miller,  M.  C. ;  "Experimental  Studies  of  Low  Viscosity, 

Immiscible  Additives  to  Reduce  the  Destabilizing  Moment 
Produced  by  a  Viscous  Liquid",  Proceedings  of  the  1990  CRDEC 
Scientific  Conference  on  Chemical  Defense  Research, 

CRDEC-SP-034 ,  August  1991. 

-A 

Weber,  D.  J. ;  "Effect  of  Surface  Roughness  on  the  Creation  of 
Viscous  Liquid-Fill  Flight  Instabilities",  Proceedings  of  the 

1990  CRPEC  Scientific  Confertnce  on  Chemical  Defense  Research, 
CRDEC-SP-034,  August  1991. 

t» 


26 


Selmi,  M.  and  Herbert,  T. ;  "Two-Fluid  Flow  in  Spinning  and 
Nutating  Cylinders" ,  Proceedings  of  the  1990  CRDEC  Scientific 
Conference  on  Chemical  Defense  Research,  CRDEC-SP-034 ,  August 
1991. 

Herbert,  T. ;  "Computer  Visualization  of  Flows",  Proceedings  of 
the  1990  CRDEC  Scientific  Conference  on  Chemical  Defense 
Research,  CRDEC-SP- 034 ,  August  1991. 

Li,  R.  and  Herbert,  T. ;  "Numerical  Study  of  Unsteady,  3D  Flows 
in  a  Spinning  and  Nutating  Cylinder",  Proceedings  of  uhe  1990 
CRDEC  Scientific  Conference  on  Chemical  Defense  Research, 
CRDEC-SP-034,  August  1991. 

Miller,  M.  C. ;  "Elimination  of  Viscous  Liquid-Fill  Flight 
Instability  by  Means  of  Lower  Viscosity,  Immiscible,  Liquid 
Additive",  AIAA  29th  Aerospace  Sciences  Meeting,  Paper  No. 
AIAA-91-0679,  January  1991. 


# 


21 


CRDEC  BASIC  RESEARCH 
IN  FLUID  DYNAMICS 


CONDUCT  EXPERIMENTAL  AND  THEORETICAL 
INVESTIGATIONS  OF  PHYSICAL  PROPERTIES,  DYNAMIC 
CHARACTERSTICS,  AND  NEW  PHENOMENA  RELATED  TO 
THE  DELIVERY  OF  CHEMICAL  PAYLOADS  ASSOCIATED 
WITH  ANTIMATERIEL,  SMOKE  &  OBSCURATION  AND 
FLAME  &  INCENDIARY  APPLICATIONS. 


FLUID  DYNAMICS 


_ TASK _ 

1.  Fluid-Filled 
Projectile  Flight 
Instabilities 


2.  Fluid  Rheology 


_ OBJECTIVE _ 

Perform  experimental  and 
theoretical  studies  to  under* 
stand,  predict  and  control  flight 
instabilities  of  liquid-filled 
projectiles  for  chemical 
munitions. 

Perform  experimental  and 
theoretical  studies  to 
establish  Improved  techniques  for 
techniques  for  determining 
rheological  properties  and 
dynamic  behavior  of  chemical 
compositions. 
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CRDEC  EMPHASIS  ON 
LIQUID-FILLED  PROJECTILE  RESEARCH 


A0332-*103-LS1  1-12327 


«*  tLm'm 


FLUID  DYNAMICS 
(FY91  ACCOMPLISHMENTS) 


TITLE:  FLUID-FILLED  PROJECTILE  FLIGHT  INSTABILITIES 

*  Theoretical:  (CRDEC  research  contract  with  The  Ohio  State 
University) 

*  Developed  advanced  computer  graphics  techniques  to 
illustrate,  Interpret  and  validate  internal  flow  - 
includes  still  and  animated  displays  (Anatomy  of 
Liquid-Filled  Projectiles). 

*  Analysis  performed  for  two  immiscible  fluids  with  large 
differences  in  viscosity  with  the  lower  viscosity 

fluid  present  in  small  amounts. 

*  Evaluated  effects  of  transient  conditions  on  creation 
of  liquid-fill  Induced  flight  instabilities  and  compared 
wnh  experimental  results. 
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FLUID  DYNAMICS 
(FY91  ACCOMPLISHMENTS) 


TITLE:  FLUID-FILLED  PROJECTILE  FLIGHT  INSTABILITIES 

*  Experimental:  (Research  studies  at  the  CRDEC) 

*  Conducted  experimental  investigations  using  CRDEC 

Laboratory  Test  Fixture  For  Non-Rigid  Payloads: 

*  Effect  of  two  immiscible  fluids  with  large  differences 
in  viscosity  on  reduction  of  viscous  liquid-fill  flight 
instabilities.  (Concept  evolved  In  collaboration  with 
Dr.  D.  Joseph,  University  of  Minnesota.) 

*  Effect  of  internal  surface  roughness  of  payload 
compartment  on  creation  of  viscous  liquid-fill  flight 
instability. 

*  Operation  Desert  Shield  program  to  rapidly  field 
artillery  projectile  having  special  liquid-fill 
similar  to  those  Investigated  previously  under  the 
CRDEC  Basic  Research  Program. 
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EFFECT  OF  ADDITIVE  VOLUME 
WITH  100K  CS  LIQUID-FILL 
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FLUID  DYNAMICS 

(FY91  ACCOMPLISHMENTS)  continued 
TITLE:  FLUID-FILLED  PROJECTILE  FLIGHT  INSTABILITIES 


*  Completed  Instrumented  flight  tests  of  artillery 
projectiles  to  validate  experimental/theoretical 
results  of  partial  fill,  viscoelastic  fill  and 
Immiscible,  low  viscosity  additives. 

*  Modified  Epicyclic  Theory  to  include  effect  of  viscous 
liquid-fill  to  predict  flight  stability  using  either 
theoretical  or  experimental  data. 

*  Co-sponsored  "Workplace  on  Problems  of  Rotating  Liquids" 
at  the  Army  High  Performance  Computing  Research  Center." 
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INSTRUMENTED  FLIGHT  TESTS 
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INSTRUMENTATION:  Yaw  Sondes  and  radar 
LOCATION:  Dugway  Proving  Ground 
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FLUID  DYNAMICS 
FUTURE  DIRECTIONS 


*  Continue  to  investigate  methods  to  reduce  or  eliminate 
viscous  liquid-fill  flight  Instabilities. 

-  Immiscible,  low  viscosity  additives. 

-  Longitudinal  baffles. 

*  Evaluate  effect  of  non-cyllndrlcal  and  eccentrically 
located  payload  compartments  on  producing  liquid-filled 
flight  instabilities. 

*  Study  neutrally  buoyant,  second  body  phenomenon  related 
to  flight  stability. 

.*  Consider  flight  stability  effects  of  novel  fluid-fills 
Including  slurries,  powders,  viscoelastic  fluids,  etc. 

*  Exploit  unusual  behavior  of  novel  fluids  to  increase 
control  of  flight  vehicles  and  other  fluid  dynamic 
devices. 
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DESPIN  VERSUS  KINEMATIC  VISCOSITY 
SPIN  FIXTURE  DATA 
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AREAS  OF  CONCERN  IDENTIFIED  AT  ROUNDTABLE 
ON  LIQUID-FILLED  SHELL,  20-21  SEPT  1984 
(not  prioritized) 

Effect  of: 

Changes  in  internal  geometry  (damp  out  inertial  modes) . 

Interaction  of  internal  flow  conditions  with  external  aerodynamics. 
Unconventional  internal  geometry  and  configurations. 

Partial  solid/Partial  liquid  payloads  including  porous  media. 
Non-Newtonian  fluids. 

Spin-up  and  cone-up  (Transient  behavior  of  fluid  and  shell's  response). 
Mixing  of  several  fluids. 

Chemical  Reactions  in  Multi-component  fluid  systems. 

Linear  vs  non-linear  analysis. 

Thermal  effects. 

Understand  and  Predict: 

Hydrodynamic  and  aerodynamic  stability  and  their  coupling. 

Theoretical  basis  of  bias  due  to  frequency. 

Sensitivity  of  shell's  response  to  cavity  aspect  ratio. 
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Work  To  Do/Gaps: 


Experimental  (Laboratry) : 

Partially-filled  container. 

Verification  of  Whiting's  experimental  results  at  low  angles. 
Pressures  at  low  Reynolds  numbers. 

Velocities. 

Despin  and  yaw  moments  at  ail  Reynolds  numbers. 

Experimental  (Plight  Test) : 

Yaw  Sonde  for  different  Reynolds  numbers. 

Identify  capabilities  and  limitations  of  current  facilities. 
Analysis: 

Establish  systematic  approach  for  all  general  fluids. 

Exploit  all  methods:  linear,  perturbation,  asymptotic,  energy 
finite  difference  and  finite  element. 

Examine  bifurcation  parameters. 

Analyze  spin-up  using  "matched  asympototoc  expansions. 

Analyze  conditions  for  Stewartson  instability  due  to  coning  match. 

Computations: 

Eigenf requencies  for  partial-fill. 

CRAY  II  for  coupling  of  aerodynamics  and  fluid  dynamics. 

Finite  difference  and  finite  element  with  adaptive  grids. 

Comparisons : 

Analysis  with  computations. 

Analysis  with  experiments. 

Develop  standardized  units  and  parameters. 

Determine  ranges  of  parameters  of  interest  to  Army. 

Large  angle  (90  degree)  coning  experiments. 


Future  Directions  - 
Liqu 1 d-Fl I  1 eo  Projectiles 


Exper imental : 

Non-cyl lndrical  containers  (all  viscosit 

Longitudinal  baffles  to  eliminate  flight 
of  viscous  liquid-fills. 

Immiscible,  low  viscosity  additive  to  el 
instabilltes  of  viscous  fills  (transient 

Effects  of  non-Newtonian  fluids. 


Theoret ica 1 : 

Purely  theoretical  method  that  can  handl 
(one  stop  approach). 

viscoelastic  fluids. 

Computat lonal : 

Launch  transient  effects. 

Liquid-Filled  Projectile  Design  Handbook 


les) . 

instabi 1 ities 

iminate  flight 
effects) . 


e  al 1  viscosities 
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Analysis  and  Visualization 
of  the  Flow  in  a 

Spinning  and  Nutating  Cylinder 


Thorwald  Herbert 
Rihua  Li 
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Department  of  Mechanical  Engineering 
The  Ohio  State  University 


Supported  by 
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Outline 


•  Introduction  -  History 

•  Analytical  Studies 

•  3D  Spectral  Navier-Stokes  Solver 

•  Eigenfunction  Expansions 

•  Computation  of  the  Moments 

•  Flight  Simulations 

•  Extensions  (Li,  Selmi) 

•  Experimental  Flow  Visualization 

•  Computer  Visualization 

•  Summary 
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FIG .  3  •  YAVJSONDE  DATA  FOR  XM76I 
STABLE  FLIGHT  (WP  SOLID) 


FIG. 4.YAWS0HDE  DATA  FOR  XM761 

UNSTABLE  FLIGHT  (WP  LIQUID) 
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Axial  Vgto C*TV 
*R«.  *  > 
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HoMgNTS 
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moment 

4o  CorloLs  forces 


UoTE  :  H  ..  =  oce/') 


X  Miller  (19B2)  c^|0. 
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Cy>STtNOi  COPE-S 
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NAviER- STOKES  SOLVER 


•  deviation  ^rovr>  solid -bodu  rotation 
Is  governed  by  a  smo.il  p am.vn el-cr 

-*■  good  approximation  j$r©m  M nest' 
equations 

-*■  $e«»3  ftsurier  modes  in  <f 

•  loco  Reynolds  number  range 

-*•  ^evo  C**»eb ijshev  modes  in  rsE 


•  Spectral  Collocation  method  : 


t  M  N 
\/<s  j  j  j 


del  Hi *i  n*t 

similar  -$or  %y  $  p 


atm«  ‘RCt'rt  Z«  C I  > 


< 


—  Symmetries  are  exploited 

—  boundary  Conditions  Implicitly 
Satisfied  by  expansion  tunetions 

“  interior  C  sine)  collocation  points 
to  avoid  corner  problem  an©| 
Spurious  pressure  term S 


SPECTRAL  COPE 
^  MoVw'oal  variables 

cz=S"  '  »  \  ^qtb'Touc  6vi«Vfcvn 

howlmearJ  13  3 

o£  d«mev>«»on  A*L.*M*N 

Cause  elitninaKon  + 
mcdlf-Ied  New+ori  iVoved'iow 

TjpicaU^  :  ALHN*^©® 

•  'Rapid  Covw/tiCQttoCe,  (NevcrtW) 

•  Pap  Id  cdnv&’oje^cc-  *  w  <f  2 

•  Semi  -  aualM+i'cal  soIu»Vi’©k> 

•  Small  data  volume 

•  Ions 

sVabllifu  awaits I5 

uweffcadvi  problems 

•  Apollo  PNSoo  (  Cravj  IS  / 2.000  ) 

S-S~s  (SOO)  9o  win 
M-3-H  (HZ)  ta3vn\r> 
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Analysis  of  uquip  moments 

4®*“  en^inte»r»wa  design  f 
•  quasi  -  steady  motion 


A^  Navier-Stofees  >  3D4  nonliVitaf 

*Re  *  1000  4  s  -  40  sec  /solut  ion 
V'erbert  e  U  ( speetml  ceole  5 


©  }  ’Boundary  layer  approximation 
“Re  Si  looo  ,  ? 

Mu-rphu 


C>  Perturbation  expansion  in 

fi>  *  §  ***  ©  (  large  c/a  *  *j  > 

all  *Rc  4  o.off  sec  /solution 
Merbert  i>  Ui 

’O)  Spatial  c^^en^unction  expansion 

Oineatieed  ,  M*CMg>  only  ) 

"Re  £  1000  >  io  -  i^oo  aec/solu+ion 
Mall  4  Sedw«y  >Gwrbc.<r 


Conclusions 

•  ,re*ult.s  o^'ree  at  *Re*'©oo 
to  lOiHim  cl  ^tw  percent 

•  UneonBation  m  e  •§«»«© 

Causes  Svnoll  error 

•  time  per  solution  »s  acceptable 

•  'time  per  simulation 

•s  considerable  C  »ooo  solutions') 

can  we  do  better  ? 
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Governing  Equations 


•  Navier-Stokes  equations 

•  Nutating  (aeroballistic)  coordinates 

•  Linearized  in  £  =  “Sin0 

co 


V  •  v  =  0 

v  +  2x  x  v  +  Vpd - ~V2v 

9(j)  yd  Re 


=  —  2r  Tre, 


where 


v2=  JL  1A  -A 

9r2  r  dr  y2  9(j)2  9z2 

Tr  =  —  —  sin0cos(J) 
co 


(^r  >  ^(h»  1  )  » 


T(h  =  —  sinGsincb 
v  co 

x  =  Acos0 
co 


7B 


•  Take  the  curl  V  x  of  the  momentum  equations 

•  Take  the  curl  V  x  of  the  resulting  equations 

3  1  9 

•  Apply  (-r - Vz)  to  the  resulting  equations 

d(j)  Re 


3(j>2 


V2v  -  ^-4-V4v  +  Av6v  +  4(1  +  tz)2-^-  = 


Re  90 


Re 


dz 


0 


•  Introduce  Fourier  series  for  v 


oo 


v=  E  *„(/-,  z)e 

n  =  -  oo 


ini j) 


•  For  the  moments,  we  need  only  vz  1?  v^q 


V2  r— 7 A  /s 

v7  +  — V  v7 
z  to  z 


1a  9 

A  ~  yi/1  ,  „  \l  4 


TVbvz  -  4(1  +  xzr — f  =  o 
2  9z2 


where 


v2  =  ~  +  --|- 

dr2  r  dr 


1  _9i 

3z2 
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Scaled  Equations  and  Boundary  Conditions 

•  Introduce  vz  =  vzl 

r  —  qr  ,  z  =  qz 

q 2  =  iRe  ,  q  =  (1  +  i)^Rel2  , 

co*\plcv>  lo/rge  ? 

a2  ia  i  32 

dr2  r  dr  r2  dz2 

•  Governing  equation  for  vz 

V2v2  -  2V4vz  +  V6vz  -  4(1  +  Tz)2^-  =  0 

dz 

s  \  xHi  -  "PPE 
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Boundary  conditions  at  the  ends: 


y 


where 


ei 


Boundary  conditions  at  the  wall: 


v 


2 


0 


2(1  +  %  )[■ 


3vz 

dr 


Vrv,  -  2- 


V  1 

dr  2 


~  ?\~2 


1  f,2  lfr4  2ft?32 

.  V]  vz  +  —  Vj  v2  +  —  Vj 

r  r  r  dz 


dr  dz 
v 


~2 


1  3v2 

+  T~1 
r  cr 

2(1  +  xz)e 
iq 


* 


where 


a2  +  _i__a_ 

a? 2  f  dr 


SPATIAL  giggN^HHCrioNS 

Assume.  vz  =  W<r>  F<«) 

V*W  -BW*0 

+hen  W<?>«r,<rBT) 

F<5)  »  £  C;  COS  a;Z 

*e  Ok  solution  of+he  Poe  «£  o>  are  roots  of 

a*  -  oJ*<3B-2^a*CSB4-*i*B  +  i-H(i  +  r^3 

-  (B*  -2s*  +  B>  »  O 

Pne  end- wall  conditions  are  satisfied  »f 
(ef-a^KSB-i-c^-a*)  a,4an  (a,q»'p  a 

4  a3  ’ I -oj-aj  )«x4a»i  < a5+eui  (a^) 
♦  (a£-  a*  >  <zB  - 1  -  c£ -  a*  >  a4+an  (os<\rj )  a,4an  fa<j  ) 

=-  o 

a«ol  +<oo  of  c •  are  proper  lu  chosen 
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Re-  1.000  Axial  Eigenfunction  n-0 

tauz-  0.200000  -  real  part 

eta-  4.000000  -  imaginary  part 
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tauz-  0.200000  -  real  part 

eta=  4.000000  -  imaginary  part 


se«>es  solution  fo*  vzi 


j 


*1*0 


—  N 


*0®.  4  ermine,  A*m  from  bounda.ru 
conditions  on  the.  s*dc  wo\l  ^so 

O  by  collocation  method 
2)  by  least  souartS 
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nJess  Momer 


Convergence 


[•ill] 


date:  11/10/88 


H'x/(pC32a5) 


Complete  Fill:  Roll  moment  vs.  Aspect  ratio 
Q  =  2,  t  =  0.083333 


up  -h>  2Zooo 


* 


90 


Conclusions 

•  Two  PT>€  vr)  is  a  s^ood  basis  for 
vviore  efficient  wiowtvi+  calcu lotion 

•  knowledge  of  tbe  analytical  structure 
of  v^i  replaces  +rial-anel-crror 
approach 

•  Spa+ial  eigen-functions  im  Z.  a. rt*oun 
alternative  -to  H  SQ  ( i»  r  ) 

•  “Results  of  eigenfunction  expansion 
®<gre<e  with  spectral  NS  reeultsat 
small  £  (  I inearisat  ion  ) 

•  Soth  methods  work  easily  at  Re*  loo* 

•  "Roots  c.<x\cu.\aJreA  up  4c  *ke  «  \oC 

•  Solution®  calculated  up  to 
*Re  *  2S  ooo 


MoHgNTS 

•  SM^aee  approach  (  ata.Holo.yol  >  .• 

JS  s4re««es  on  instate,  wall 
S 

pressure 

Veloc^  gradients 
c*'H*’'rftv,c *>  of  lav«e  numbers 

•  Volume  approach  : 

ill  Coriolis  4-e.rms  in+Ht  fluiol 

y 

Tlosenblats  relations 
av«al  z>  aii  m  u  fh  a  I  velocities 

tbr  given  fields  »44te  volume  approach 
provides  more  accurate  moments 

U i+h  +hc  volume  approach  4  ee4i‘mates 
for  4-he  moments  can  be  obtain eol 
from  +lte  0unalvj4ica,l  res u  14s 
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CALCULATION  OF  MOMENTS 


i 

*  ^  OO*  ©  J  J  |  V L'T^COS’f  dkdfdz 

•1)0  0  t  V 

H*R» 

*  !y«;"©  J  1  J  v„  r‘drelt <**■ 

|  '  O  0  a. 

&  V|  llv  • 

+2  geos©  $  J  J  vz+'s;„n  dt^dx 


•loo 


*  M„  +av>  © 


•  «*P 


Le*  21  vc^e. 

4+»ew  we  need!  onl^ 


v*i<r,,z'>  v<f0<’,*.*> 

'  -  i 

<s*\a\l  e^ter 
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log  (Re) 
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fui<RHT  INSr*SlUTt£S  : 

•  *««.4oiiancft  m.H>  inerKaf  waves 

C  \*:§h  *^*3  v>o  l©l« 

Kwmbtr  ) 

•  Vtceoue  secondary  §too0 

C  vweoU  uvn  *RcAjiAo((»l-9  Muwibfcf') 

Flight  -Simulation  ; 

Va-w^Wn^  Wolfe  10b«,'r)*a»wjoC 

shell  dabcw 

ae'rooitjwawie  da+a.  C+able) 

C|oo^j load  da>a 

It^uxol  vwomewfa  C+abie) 
lo/wndx  dala 

» n+t^trd’ion  3o - 1005  ,  M«  2- io‘?5 


so sr  t 


FLi<5HT  SiMULATlGN  ••  Avi^lt  of  O.M<K^ 


9  0S 


gy  ^gRtHgKTT^L  STUPtgS 

Goal  Show  •feasibUi'bj  of -flow 

Visual  iWt'ons  et  pfo^criu 
Scot.)  col  4*oolcl 

T?tc+fid’*on*  :  low  CO«t*  (-4  S'oo  ) 

use  ava>  fable  t^ulamcnt 
uhlitt  ''cheap*  lowboy 

^  Stniof  S+moIcwI’  Pyoleeh : 

'Pay*©!  'Pi’crponl- 
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Mean 

NSY_5  /  NSY_i  Velocity  Field  -  Type  o 


Re» 

Theta* 

Tail* 

Lambda* 


20.000 
20.000 
0.1667 
4 . 3680 


Cut  at  Phi-  0 

Scale:  - -  0.00200 

Points:  6  radial 

5  azimuthal 
S  axial 


* 


h  I 


NSY_5  /  NSY_1  Velocity  Field  -  Type  0 

Re-  300.000  Cut  at  Phi-  0 

Theta-  20.000  Scale:  - -  0.05000 

Tau-  0.1667  Points:  6  radial 

Lambda-  4.3680  5  azimuthal 

8  axial 


NSY_5  /  NSY_1  Pressure  Field  -  Type  o 

Re-  300,000  Cut  at  Phi-  0 

Theta-  20.000  Levels;  0,00500 

Tau-  0.1667  Points:  6  radial 

Lambda-  4.3680  5  azimuthal 

8  axial 


NSY _5  /  NSY.J 


T©4*| 

Pressure  Field  -  Type  o 


Re-  300.000 
Theta-  20.000 
Tau-  0.1667 
Lambda-  4.3680 


Cut  at  Phi-  0 
Levels:  0.04000 
Points:  6  redial 

5  azimuthal 
8  axial 


116 


CONTOUR  LEUELS 


i  “  «■*»  • 
*.  >_  L1  w 


-II  .  t]  1  (II  ill 

-o.  nor.no 
0.00000 
o.onsoo 
cr,  mono 

O.OlbOU 

n.  n?pnn 

o, n?r.nn 
0,00000 
O.OObOO 


t 


♦ 


Bl  ank 


W 


120 


ROTATING  FLUIDS  WORKSHOP 
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University  of  Minnesota 


i  Q 


QP  =  <*(i  -  i)Q 

+aUs(H  +  b)2nR 
If  Qp  =  aQ, 

Us  =  'jQ/(H  +b)2nR 


0  ~  i)Q 


SMALL  PARTICLES 

aj 

—  =  a 
a  2 

%  OUT  OF  SLOT 


iQ 


=  1  -  7(1  -  a2) 


%  REMAINING  =  7(1  -  a2) 
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V  •  q  =  0 


— q  +  -Roq  •  Vq 

+2 k  x  q  =  -VP  +  EV* q 


e  = 


Pp  ~  Pf 
-  • 


2a2n 

9  v 


B.C.  q  —  Uvy  —  Uw ft 
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VERTICAL  B.  L.  SCALE  ~  {u/UH2)l^H 

~  El/2H 

ENTRY  LENGTH  ~  SQ/Au 

Ro  ~  E^2  N.  L.  EFFECTS 
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Moment  Exerted  by  a  Visco 
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Ballistic  Research  Laboratory 


LINEAR  LIQUID  MOMENT 


s 

CD 

0 

CL 

03 

r-> 

+ 

03 

03 

5 

o 

C 

5 

_i 

“O 

03 

o 

c 

0 

03 

CC 

4- 

3 

0 

LSM 

_j 

o 

o 

L__l 

CD 

o 

OJ 

-e- 

CM 

o3 

1 

0 

c 

o 

03 

E 

o 

c 

Z3 

4 


133 


Classical  Method 
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proper  selection  of  ck,  Xk,  ai 
the  b.c.  of  the  R,.  functions 


LIQUID-FILLED  CrLINMR  THEOty 
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H64  HflLL-iEDMEy-GERBER  (LNS) 


MOMENTUM  DERIVATIVE 
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for  A  =  4.5,  t  =  0.1.  Comparison  with  MHUr  data. 
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Partial  Pill/Central  Rod 
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0  P  0 
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Partial  Fill:  CLSM  vs. 
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Figure  8.  Ratio  of 
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NUMERICAL  STUDY  OF  UNSTEADY  3D  FLOWS 
IN  A  SPINNING  AND  NUTATING  CYLINDER 


* 


R.  LI  &  TH.  HERBERT 


RESEARCH  FUNDED  BY  CRDEC 

DEPARTMENT  OF  MECHANICAL  ENGINEERING 
THE  OHIO  STATE  UNIVERSITY 
COLUMBUS,  OHIO 

APRIL  1991 
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V 


PREVIOUS  WORK:  STEADY  STATE  ONLY 

*  FDM  (USE  ARTIFICIAL  COMPRESSIBILITY) 

VAUGH,  OBERKAMPF  &  WOLFE  (1985,  JFM) 

*  FDM  IN  R  &  Z,  FOURIER  METHOD  IN  <j>. 

STRIKWERDA  &  NAGEL  (1988,  J.  COMPUT.  PHYS) 

*  SCM  (SPECTRAL  COLLOCATION  METHOD) 

HERBERT  &  LI  (1990,  J.  AIAA) 

PRESENT  WORK:  UNSTEADY  (SPIN-DOWN) 

Q  AND  0  ARE  FIXED,  WHILE  co(0  IS  GIVEN. 

n 


G) 


Fig.  1  A  sketch  of  a  nutating  and  spinning  cylinder.  The  frame  ( x-y-z  )  is  rotating 
with  respect  to  the  fixed  frame  (X-Y-Z)  with  the  angular  velocity  Q.. 
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VELOCITY- VORTICITY  FORMULATION 

-|-0  +  (V  V)  O  =  (O-V)V  +  vV20  +  20  VV  -  20  (1) 

ot 

V2V  =  Vx  O  (2) 

v  :  KINETIC  VISCOSITY 
Q  :  ANGULAR  VELOCITY  OF  NUTATION 
V  :  VELOCITY 

O  :  VORTICITY  (O  =  -VxV) 

REFERENCE  QUANTITIES  FOR  NONDIMENSIONALIZATION: 
RADIUS  a,  DENSITY  p,  AND  <oref 

NUMERICAL  METHOD 

*  FDM  IN  R  &  Z  (UNIFORM  GRID:  20  PTS  IN  R;  40  PTS  IN 
Z); 

*  FOURIER  METHOD  IN  <J>  (6  FOURIER  MODES); 

*  VISCOUS  TERMS  -  MULTISTEP  IMPLICIT; 

*  CONVECTIVE  TERMS  -  MULTISTEP  EXPLICIT; 

*  2ND  ORDER  ACCURACY  IN  SPACE  &  TIME; 

*  ITERATION;  LINE  GAUSS-SEIDEL  RELAXATION 

WITH  ZEBRA  PATTERN  (VECTORIZATION). 
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CALCULATION  OF  MOMENTS  : 

VOLUME  INTEGRAL 
SURFACE  INTEGRAL 

VOLUME  INTEGRAL  APPROACH 

M  =  f  r  x  2^-dV  =  vJ  +  AxJ  (3) 

J  Dt  3t 


IF  WE  DECOMPOSE  THE  VELOCITY  FIELD  INTO 

V  =  Vr  +  v*  (4) 

Vr  :  RIGID  BODY  MOTION;  Vd  :  DEVIATION. 

THEN 


j  =  y  +  3d 

(5) 

WHERE 

=  (A  23  ,  -A  f  3  ,  A  d 2  ) 

(6) 

INTEGRALS  (DIMENSIONLESS) 

A  f3  =  jvfr  cos<prdrd  <pdz 

(7) 

A  23  =  jVfr  sinfyrdrd  §dz 

(8) 

A  12  =  ~jV$r2drd  §dz 

(9) 

(INDEED,  ONLY  NEED  FUNDAMENTAL  OF  Vd 
TION  OF  THE  MEAN  OF  Vd). 

AND  DISTOR- 

f 
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DIMENSIONLESS  COMPONENTS  OF  MOMENT 


Mx=2Adn 

+  2eA  f3  /  tan0 

(10) 

My  =  -  2/if 3 

+  2e(A  12  +  A  23  /tan0) 

(11) 

Mz  =  2A  f2  +  r|7tcb(r ) 

+  2£Af3 

(12) 

« 

—  o  2  ?  1 

My  =  My  +  £C0(/)rj7t  +  e^r|7c( — rf  +  — )/tan0 

3  2 

03) 

PARAMETERS:  e  =  Qsin0/corf/ ,  r|  =c/a  (ASPECT  RATIO) 


WHERE 

(1)  COMPONENTS  DUE  TO  THE  UNSTEADINESS  OF  RELA- 
TIVE  MOTION  (THE  TERM 

RIGID  BODY:  (  0,  0,  T| 

DEVIATION:  (  TA  f3 ,  -  2 A  f 3 ,  1A  f2  ) 

(2)  COMPONENTS  DUE  TO  THE  CORIOLIS  ACCELERATION 

RIGID  BODY:  (  0,  ec o(r  )t|n,  0  ) 

DEVIATION:  (  2eAf3/ tan0,  2£(Af2  +  Af3/tan0),  2sAdn) 

(3)  COMPONENTS  DUE  TO  THE  CENTRIPETAL  ACCELERA¬ 
TION 

RIGID  BODY:  (  0,  e2t|n( — 2_r|2  +  -)/tan0,  0  ) 

3  2 
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moment 


RESULTS  AND  DISCUSSIONS 

(1)  CALCULATION  OF  STEADY  STATES 

*  START  FROM  RIGID  BODY  MOTION 

*  COMPARE  WITH  SCM  (5  FOURIER  MODES  in  <j), 

5  CHEBYSHEV  POLYNOMIALS  IN  RADIUS, 

5  CHED  \SHEV  POLYNOMIALS  IN  HALF  LENGTH.) 

*  Q  =  500  RPM,  co  =  3000  RPM,  6  =  20  DEG,  <\ef  =  co, 

rj  =  c la  =  4.368,  Re  =  a2co/v. 

*  0.3  TO  2  MINUTES  IN  CRAY  Y-MP  FOR  Re  <  200. 
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(2)  SPIN-DOWN  OF  60  KCS  FLUID 

*  Q  =  500  RPM,  0  =  20  DEG,  Tj  =  4.4955. 
core/  =  20  =  1000  RPM. 

"INSTANTANEOUS  REYNOLDS  NO."  :  Re  =  a2co(f)/v. 


8000 

7000 

^  6000 

a 

£;  5000 

"  4000 
(0 

c  3000 
a 

M  2000 
1000 
0 

0  5  10  15  20  25  30  35 

time  (sec) 


(A)  CALCULATE  STEADY  STATES  AT  DIFFERENT  SPIN 
RATE. 

(B)  START  FROM  STEADY  STATE  AT  t  =  0,  USE  THE 
GIVEN  co(0  AS  TIME-DEPENDENT  BOUNDARY,  RUN  THE 
"QUASI-STEADY"  SPIN-DOWN. 

*  TIME  t  =  0  TO  t  =  3000 
At  =  0.02,  150,000  TIMESTEP,  5  HRS.  IN  CRAY. 


Time  History  for  60KCS  Fluid 
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Amplitudes 


EXAMPLE  FOR  STEP  (A):  STEADY  STATE  AT  (.)  =  6800  RPM 
A  ,  :  FUNDAMENTAL  OF  Vf(r  =  0.5,  z  =  0,  0  =  0) 

A0  :  DISTORTION  OF  MEAN  FLOW  OF  V^(r  =  0.5,  z  =  0,  <]>  =  0) 


0.0  1.0  2.0  3.0  4.0 

Time 


FLUID  DEVELOPS  FROM  RIGID  BODY  MOTION  TO  STEADY 
STATE  RAPIDLY. 


integrals 


THE  INTEGRALS  OF  UNSTEADY  ("QUASI-STEADY")  ARE 
THE  SAME  AS  THOSE  OF  STEADY  STATE. 


CONCLUSION  :  THE  UNSTEADY  FLOW  FIELD  ADAPTS  VERY 
QUICKLY  TO  THE  SLOWLY  CHANGING  BOUNDARY  CONDI¬ 
TIONS. 


0  2000  4000  6000  8000 

spinning  rate  (rpm) 
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Mx  and  My  OF  UNSTEADY  ("QUASI-STEADY")  ARE  THE 
SAME  AS  THOSE  OF  STEADY  STATE. 

Mz  OF  UNSTEADY  ("QUASI-STEADY")  IS  BELOW  THAT  OF 
STEADY  STATE. 


Steady  and  Unsteady  Cases  for  60KCS  Fluid 


THE  CONTRIBUTION  OF  TIME  DERIVATIVES  OF  INTEGRALS 
OF  UNSTEADY  ("QUASI-STEADY")  TO  THE  MOMENT  COM¬ 
PONENTS  ARE  NEGLIGIBLE. 

THE  CHANGE  OF  SPIN  RATE,  6(0,  HAS  CERTAIN  EFFECT 
ON  THE  RESULTS  OF  Mz  (NOTE:  T|ti6(0  REPRESENTS  THE 
UNSTEADY  RIGID  BODY  MOTION). 
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SUMMARY 


(1)  FOR  MOMENT  M  =  -§-J  +  Cl  x  J,  THE  VISCOUS  FLOW 

ot 

P) 

HAS  LITTLE  EFFECTS  TO  THE  TERM  ~J,  BUT  IMPORTANT 

9 1 

CONTRIBUTION  TO  THE  TERM  RELATED  TO  CORIOLIS 
ACCELERATION  WHICH  IS  INCLUDED  IN  THE  TERM  QxJ. 

(2)  THE  MOMENTS  OBTAINED  FROM  THE  <  "QUASI- 
STEADY"  SPIN-DOWN  EXPERIMENTS  CAN  BE  REGARDED 
AS  THOSE  FOR  STEADY  STATES  IF  THEY  ARE  CORRECTED 
BY  TAKING  THE  UNSTEADY  RIGID  BODY  MOTION  (THE 
TERM  TiTub(O)  INTO  ACCOUNT. 

(3)  FOR  STEADY  FLOW,  Mz  =  -Mx  sinG  +  M2cos0  =  0 
BUT  FOR  UNSTEADY  FLOW, 

Mz  =  -AffsinG  +  M* cos9  =  -2(A  23  sin0  -  Adl2  cos0)  +  \n( o(t  )cos0 

(4)  THE  FORMULA  OF  MOMENT  CAN  BE  EASILY 
EXTENDED  TO  THE  UNSTEADY  CASE  THAT  BOTH  Q.  AND  co 
ARE  CHANGING  WHILE  0  IS  FIXED  (FOR  EXAMPLE,  "SPIN- 
UP"). 


Modelling  Non- Newtonian  Beliavi'or 
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Obvious  non  —  (Nitwfoni'an  jbejtat/ior  •for  krjh  viscosity 
Silicone  200  :flu*d s  6for  example ,  V'?6okcs) 

S.  Rosenblat  et  aj.  OW) 

Analytical  Work  (.  perturbation  method) 

Numerical  Stud  res  CFEM) 

No  measured  data  from  experiments  available. 
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P-  P.  T/tus  (118*1,0*0 EO 

measured  Silicone  2.00  -fluids 
Used  the  -formal#,  of-  #t»ffe.renti'a(  model  derived  by 
Rosen bUt  to  calculate  tfu  clesp»n  moment  —  good 
agreement  coes  not  achieved. 


OiHiri nti»l  r»od»1 


I  2  3 

LOO  R# 


Problem  :  tun?  terms  are  missing  in  J?ose*bU+'s 
7$t>r  mala . 

The  Corrected  flo rmuta.  seems  better  for  loo  Ires. 

Future  work-  1.  Check  other  Ugh  vtecositfes; 

2.  3-D  rucmenW  studies . 
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Non- Cylindrical  Containers 

Previous  Analytical  and  Experimental  Work 
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Generalized  geometry,  Zonal  mesh 
Volume  integral  method  for  Cxrm^lsm  added 
Turbulence  and  separated  flow  models 
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Computational  Grid 

Truncated  Ellipsoid 
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•  Container  axis  at  interface  of  cylinder  and  endcaps  (r  =  0,  z  —  \c/a\ ) 
assigned  equal  values  (obtained  from  interpolation) 


UWISC/BRL  -  Numerical  Method 
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where  m  is  the  iteration  index  and  (3  is  a  parameter 
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Cylinder  Aspect  Ratio,  c/a 


UWISC/BRL  —  Conclusions 
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Id  matrix  accounts  for  moving  grid  zones 


USA-IN3  —  Numerical  Method 
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Implicit  time  discretization  up  to  third-order  accurate 

•  Sub-iteration  guarantees  divergence-free  solution  at  each  At 

•  ADI  and  various  Gauss-Seidel  methods  in  either  full-block 
or  diagonalized  formulation 


USA-IN3  -  Boundary  Conditions 
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USA-IN3  Result 

Vortex  Shedding  Behind  a  Circular 
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USA-IN3  Result 

ortex  Shedding  Behind  a  Circular  Cylinder  (Re#  =  2000) 
Strouhal  Number  —  .21  (experimental),  .23  (computed) 
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USA-IN3  Result 

Vortex  Shedding  Behind  a  Sphere  (Rep  =  2000) 


USA-1N3  Computational  Fluid  Dynamics  Code 
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(CRAY-2  Computer) 
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GOVERNING  EQUATIONS 
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NON-DIMENSIONAL  PARAMETERS 
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FLOW  DECOMPOSITION 
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Complete  Fill:  Roll  moment  vs.  Reynolds  Number 
77  =  4.368,  6  =  20.  t  =  0.16667 
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Complete  Fill:  Roll  moment  vs.  Aspect  ratio 
6  =  2,  r  =  0.083333 
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207 


Two  Fluids:  Roll  moment  ra 
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Two  Fluids:  Roll  moment  vs.  Fill  ratio 
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Two  Fluids:  Roll  moment  ratio  vs.  Fill  ratio 
77  =  4.5,  6  =  1,  r  =  0.008 
Re0  =  25,  Re^lO*,  p0/p ,  =  0.98 
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Complete  Fill:  Roll  moment  vs.  Reynolds  Number 
77  =  4.5,  6  —  20,  t  =  0.16667 
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Mz/(pow2a5) 


Two  Fluids:  Roll  moment  vs.  Inner-fluid  Reynolds  Number 
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Pq/i P\  ~  0*98 


214 


Partial  Fill: 


Roll  moment  vs.  Fill  radius 
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Central  Rod:  Roll 
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Central  Rod:  Roll  moment  vs.  Fill  radius 

77  =  4.5,  6  =  20,  t  =  0.08674 
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Roll  moment  vs.  Fill  radius 

Re  =  1 04,  77  =  2,  6  =  2,  r  =  0.1111 
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Roll  moment  vs.  Fill  radius 


INVISCID  ANALYSIS  (  PARTIAL  FILL) 
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CENTRAL  ROD:  Critical  fill  radius  vs.  aspect  ratio 

Coning  frequency:  r=0. 08674 
Nutation  angle:  0  =  20. 


PARTIAL  FILL: 


Critical  fill  radius  vs.  aspect  ratio 
Coning  frequency:  r=0. 08674 
Nutation  angle:  8=20. 
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PARTIAL  FILL:  Interface  shape  (r0= 0.2, 0.8,  0=0) 
Re  =  10000,  t  =  0.08333,  tj  =  3,  9  =  20 
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TWO  FLUIDS:  Interface  shape  (r0  =  0.66,  <p  =  90) 
Re0  =  30,  Re1  =  10000,  t  =  0.1,  77  =  4.5,  0  =  20 
Po/p\  =  0-001 
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Liquid  Moment  Magnitude  History 
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Experiment  Selection  Criteria 

Hardware  -  c/a,  lx  &  ly,  moment  sensitivity,  offset 

single-mode  motion,  spin  rate,  coning  rate  and  amplitude 
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Reynolds  Number 


240 


Liquid 


* 


SO  ^  O 

co  ^  a 

Os  IN  ^  IN 

o  co  oo  cd 


00  oo 

CO  CO  00 

oo  co  oo 


O  O  O  o 
•  ■  •  , 
T-H  rH  U-) 


LO  O  LO  o 

Tt<  ON  CO  ON 


241 


G 

o 

£ 

o 

s 

0; 

"0 

CO 

T3 

•  H 

3 

cr 

•  m 

hJ 

• 


242 


Time  (deg) 


Frequency  response  plot  of  phase  for  liquid  filled  cavity. 
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0.00  0.05  0.10  0.15 

Non-Dimensional  Coning  Rate  -  X 


Re=12,624  c/a=3.126  f=77.03%  h/a=0.0868  <*1.5  deg 

Improved  Stewartson-Wedemeyer  (volume  method) 
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nt  balance  sensitivity,  circuitry,  measurement, 
nic  frequency  response  of  moment  balance. 


Conclusions 
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Army  High  Performance  Computing  Research  Center 

University  of  Minnesota 


PROJECTILE,  1 55MM  SMOKE,  WP,  M825 
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FIGURE  2.  Internal  Arrangement  of  M825  WP  Smoke,  155mm  Artillery  Projectile 


CANISTEF; 
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M825A1  Payloads  Tested 


A  total  of  12  different  M825A1  payload  canisters  were  evaluated  during 
this  study  as  listed  below: 


CONFIG.  DESCRIPTION  CRDEC  EMPTY  FILLED  SIMULANT  I  (EMPTY)*  I  (FILLED)** 


NO. 

No. 

Wt.* 

lbs. 

Wt.** 

lbs. 

Wt. 

lbs. 

Slug-Ft. 

Slug-Ft. 

1 

STD 

L 

A 

#6 

23.88 

37.01 

13.13 

.055000 

.064838 

2 

STD 

H 

A 

#38 

24.26 

37.33 

13.07 

.055317 

.065125 

3 

STD 

L 

B 

#45 

24 . 06 

37.17 

13.11 

.055128 

.064895 

4 

STD 

H 

B 

#51 

24.16 

37.23 

13.07 

.055200 

. 065000 

5 

FAT 

L 

A 

#18 

24.24 

37.36 

13.12 

.055211 

.065046 

6 

FAT 

H 

A 

#?2 

24.32 

37.37 

13.05 

.055291 

.064968 

7 

FAT 

L 

B 

#5/ 

23.83 

36.93 

13.10 

.054930 

.064753 

8 

FAT 

H 

B 

#78 

24.10 

37 . 12 

13.02 

.055139 

.064858 

9 

FAT 

P 

A 

SHEET  #1 

24.81 

37.35 

12.54 

.055863 

.065095 

10 

FAT 

P 

A 

SHEET  #2 

24.83 

37.30 

12.47 

.055846 

.065089 

11 

FAT 

P 

A 

ROLL  #1 

24.56 

37.23 

12.67 

.055683 

.065133 

12 

FAT 

P 

A 

ROLL  #2 

24.41 

37.23 

12.82 

.055540 

.065030 

Where: 

STD  -  Standard  radius  felt  wedges,  Rfc2.46" 

FAT  -  Enlarged  radius  felt  wedges,  R«2.56" 

L  ■  Light  weight  felt  wedges,  Thickness  0 , 725"-0 . 780" 

H  -  Heavy  weight  felt  wedges,  Thickness  0.781"-0.850w 

A  -  Manufacturer  A 

B  •  Manufacturer  B 

P  -  Production  felt  wedges 

A  SHEET  ”  Manufacturer  A  felt  supplied  in  sheets 
A  ROLL  -  Manufacturer  A  felt  supplied  in  rolls 

*  Canister  contains  all  payload  components  except  WP  simulant. 

**  Canister  contains  all  payload  components  including  liquid  WP  simulant. 
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FELT  SHEET:  CLOTH,  FELT,  WOOL, 


DESIGNATION  W 
SPEC  2.44 
ENLARGED  2.54 

MOTE:  ALL  ENMCMSJONS  AME  M  INCHES 


PRESSED,  TYPE  l 
CLASS  12*3,  .75  THICK 
SPEC  CP-204 


TABLE  1 


1 

Properties  of  Fluorinert 

Fluorinert* 

White  Phosphorus 

1 

(25*  C) 

(43.3*  C) 

9  Density 

1.73  gm/ml 

1.73  gm/ml 

J  Viscosity 

1.50  CS 

1.50  CS 

|  Surface  Tension 

15.4  dynes/cm 

71.6  dynes/cm 

NOTE:  Blend  of  Fluorinert  FC40  and  FC72,  1  Part  FC40  to 

0.129  Part  0.129  Part  FC72  by  volume. 
Fluorinert  is  an  electronic  fluid  manufactured  by: 
Commercial  Chemical  Division/3M,  223-65-04  3M  Center, 
St.  Paul,  MN  55144-1000.' 
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SPIN  DECAY  (Raw  Data) 
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TIME  ~  SECONDS 


PAYLOAD  MOMENTS  vs  CANISTER  SPIN  RATE 


CANISTER  SPIN  RATE  (u)  ~  RPM  X  1000 


PAYLOAD  MOMENTS  vs  CANISTER  SPIN  RATE 


CANISTER 


PAYLOAD  MOMENTS  vs  CANISTER  SPIN  RATE 
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TABLE  2 


PAYLOAD  MOMENT  (MP)  -  FT-LBS 


PAYLOAD  MOMENTS  vs  CONING  ANGLE 


COMING  ANGLE  (0)  ~  DEG 
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PAYLOAD  MOMENTS  vs  CANISTER  SPIN  RATE 
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CANISTER  SPIN  RATE  (w)  ~  RPM  X  1000 


PAYLOAD  MOMENT  (Mpl  ~  FT  — LBS 


PAYLOAD  MOMENTS  vs  CONING  ANGLE 


CONING  ANGLE  (0) 


Conclusion: 


This  was  the  latest  in  an  ongoing  test  series  supporting  the  M825  Smoke 
Round  Program.  The  Frojectile  Flight  Simulator  or  "Spirr  Fixture”  has 
been  used  in  the  development  of  the  M825  Smoke  Round  and  in  evaluating 
sample  units  from  production  runs  for  consistency  in  materials  and 
manufacturing  process.  It  has  been  a  valuable  tool  for  both  practical 
and  research  programs. 
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Radial  compaction  is  more  efficient 
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NON-DIMENSIONAL  CONING  FREQUENCY 


COMPARISON  OF  AXIAL  AND  RADIAL  COMPACTION 


COMPACTION  RATIO 
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Comparison  of  BRL-Free  Gyro  and  CRDEC-Test  Fixture 

Felt  Experiment 

Compaction  Ratio:  Axial  =  1.0 
0051  Radial  =  1.04 
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Army  High  Performance  Computing  Research  Center 

University  of  Minnesota 


INTRODUCTION 


•  Highly  viscous  liquid  fills  (e.g.  100K  CS)  can 
cause  flight  instabilities. 

•  Small  amounts  of  low  viscosity  fluids  (e.g.  5% 
water)  added  to  the  high  viscosity  liquid  can 
eliminate  instabilities. 

•  Drawbacks  to  additive  appi.  ach: 

•  Two  fluids  must  be  immiscible. 

•  Low  viscosity  additives  must  have  a 
greater  density  than  high  viscosity  fluid. 

•  Smooth  canister  wall  may  reduce  surface 
shear  stresses. 

•  Achieve  results  similar  to  additive 
approach. 

•  Eliminates  drawbacks  of  additive 
approach. 
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CANISTER  CONFIGURATIONS  TESTED 

•  Standard  canister  c/a  =  4.5, 

(roughness  =  1  25  to  250  micro  inch  finish) 

•  Smooth  canister  c/a  =  4.5, 

(roughness  ~  9  to  17  micro  inch  finish) 

•  Smooth  canister  coated  with  silicone  based 
marine  wax. 

•  Smooth  canister  coated  with  teflon  based 
marine  wax. 

•  Rough  wall  canister  c/a  =  4.57,  standard 
canister  painted  with  grit  filled  safety  paint, 
(roughness  =  50,000  micro  inch  finish) 
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TEST  CONDITIONS 

Fluids 

•  1 0,000  CS  Silicone 

•  100,000  CS  Silicone 

•  10,000  CS  SMicone  +  5%  water  (Standard 
canister  only) 

•  100,000  CS  Silicone  +  5%  water  (Standard 
canister  only) 

Coning  Rates 

•  300  rpm 

•  400  rpm 

•  500  rpm 
Coning  Angle 

•  20  ° 
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TEST  PROCEDURE 

For  a  specific  coning  angle  and  coning  rate: 

•  The  canister  is  spun  up  to  =  10,000  rpm  as 
the  coning  rate  is  established  at  the 
predetermined  rate. 

•  With  the  coning  rate  held  constant,  the 
canister  is  allowed  to  spin  down. 

•  Canister  despin  is  recorded  versus  time. 

•  From  the  despin  data,  the  Total  Moment  (Mj) 
is  calculated,  M-p  =  I  cu,  where  I  is  the  axial 
moment  of  inertia  of  the  empty  canister. 

•  The  liquid  moment  is  found  from 

Ml  =  Mj  -  Mp,  where  Mp,  friction  moment, 
is  determined  using  the  same  procedure 
described  above  but  with  no  fluid  in  the 
canister,  Ml  =  0. 
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CONCLUSIONS 


•  For  typical  artillery  shell  flight  conditions  (Q  =  500  rpm,  uj  = 
6000  rpm)  the  additive  approach  yielded  the  best  results, 
reducing  the  despin  moment  by  ~  50%  (ref.  Standard 
Canister). 

•  10,000  cs  Silicone  Fluid 

For  typical  artillery  shell  flight  conditions  the  canister  wall 
roughness  had  little  or  no  effect  on  the  despin  moment. 

•  100,000  cs  Silicone  Fluid 

For  typical  artillery  shell  flight  conditions  the  despin  moments 
(ref.  Standard  Canister)  for  the  smooth  and  rough  canister 
configurations  were  decreased  ~  10%. 

•  Log(Re')  versus  Liquid  Payload  Coefficient 

Re'  >10:  5%  water  produced  smallest  coefficients,  standard 
and  rough  canister  produced  the  largest  coefficients  and  the 
smooth  canisters  produced  results  in-between. 

•  Above  6000  rpm,  the  rough  and  smooth  teflon  waxed 
canisters  produced  a  significant  decrease  in  despin  moment. 
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INSTRUMENTED  FLIGHT  TESTS  || 
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LOCATION:  Dugway  Proving  Ground 


DESTABILIZING  MOMENT 
DUE  TO  PARTIAL-FILL 
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VISCOELASTIC  PROPERTIES  FOR 
NEWTONIAN  FLUID  (DEM  WITH  12.1%  K-125) 
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SHEAR  RATE  (y)  ^  SEC 
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SPIN  RATE  (o>)  ^  i()3  RPM 
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DAN  JOSEPH’S  EXPERIMENT 

UNIVERSITY  OF  MINNESOTA,  1 985 
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NOTE:  1)  EQUAL  DENSITIES 

2)  IMMISCIBLE 

3)  REDUCED  TORQUE 


LIQUID-ROLLING  MOMENT  FOR  100K  CS 
LIQUID-FILL  WITH  WATER  ADDITIVE 
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EFFECT  OF  VARIOUS  ADDITIVES 
WITH  100K  CS  LIQUID-FILL 
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500  RPM 
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D-ROLLING  MOMENT  COEFFICIENT 
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PRISON  OF  NEWTONIAN  LIQUID-FILLS 
ADDITIVE  AND  VISCOELASTIC  FLUID 
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SUMMARY  OF  ADDITIVE  FLIGHT  TEST  RESULTS 


Fill  Additive 

Initial  Yaw 

Flight  Stability 

100K  CS 

NO 

Large 

Catastrophic  @  8  Sec 

N 

Yes 

tt 

"  8  12  Sec 

35K  CS 

No 

N 

N  8  13  Sec 

H 

Yes 

n 

"  8  22  Sec 

10K  CS 

No 

ft 

"  8  29  Sec 

it 

Yes 

vt 

Stable  Flight 

100K  CS 

NO 

Small 

Catastrophic  8  28  Sec 

N 

Yes 

I* 

Stable  Flight  (Limit  at  end) 

35K  CS 

No 

n 

Stable  Flight  (Limit  at  end) 

M 

Yes 

n 

Stable  Flight 

10K  CS 

No 

N 

Stable  Flight  (Limit  at  end) 

n 

Yes 

N 

Stable  Flight 

100K  CS 
(50%  Full) 

No 

Large 

Catastrophic  8  8  Sec 

If 

No 

Small 

Catastrophic  8  21  sec 

100K  CS 
(Viscoelastic) 

NO 

Large 

Stable  Flight 

If 

NO 

Small 

stable  Flight 

Note:  Zone  4  Charge,  850  nils  (47.5  deg.)  QE,  5000  ft  ground  elevation, 
M825  projectiles  with  standard,  two  piece  base,  large  yaw  -  10-12 
deg. ,  small  yaw  **  2  deg. 
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CONCLUSIONS 


*  GOOD  AGREEMENT  BETWEEN  THEORY,  LABORATORY  EXPERIMENTS 
AND  FLIGHT  TESTS  FOR  VISCOELASTIC  LIQUID-FILL. 

*  GOOD  AGREEMENT  BETWEEN  THEORY,  LABORATORY  EXPERIMENTS 

AND  FLIGHT  TESTS  FOR  PARTIAL-FILL  CONDITION.  ^ 

*  LOW  VISCOSITY,  IMMISCIBLE  LIQUID  ADDITIVE  ALWAYS 
REDUCES  DESTABILIZING  LIQUID  MOMENT. 

*  ADDITIVE  PRODUCED  STABLE  FLIGHTS  FOR  ALL  LIQUID-FILL  * 

VISCOSITIES  TESTED  UNDER  NORMAL  LAUNCH  CONDITIONS. 

*  ADDITIVE  PRODUCED  STABLE  FLIGHTS  FOR  ONLY  THE  LOWEST 
LIQUID-FILL  VISCOSITY  TESTED  UNDER  INDUCED  YAW 
LAUNCH  CONDITIONS. 

*  POSSIBLE  REASON  FOR  DIFFERENCES  BETWEEN  LABORATORY 
PREDICTIONS  AND  FLIGHT  FOR  ADDITIVE. 

-  Transient  effects  during  firing  may  have  delayed 
distribution  of  additive  to  side  wall  of  container. 

-  May  require  increased  amount  of  additive  (io%-15%) 
rather  than  5%  used  in  these  tests. 


\ 
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Summary  of  New  Directions  for 
Liquid-Filled  Projectile  Studies 


♦NOTE:  The  following  items  were  suggested  and  agreed  upon  by 
those  in  attendance  at  the  Workshop.  They  have  not  been  put 
under  any  specific  category  nor  have  they  been  prioritized. 


(_  Coordinate  efforts  between  ARO,  BRL  and  CRDEC  concerning  data, 
terms  and  comparison  of  results. 

Concerned  about  achieving  better  agreement  between  experimental 
^  and  theoretical  results. 

Investigate  heavier  immiscible,  low  viscosity  additives  (i.e., 
salt  water)  to  improve  transient  distribution  for  reducing 
viscous  liquid-fill  instabilities. 

Explain  linear  and  non-linear  effects. 

Conduct  bifurcation  studies. 

Determine  density  of  felt  wedges  when  saturated  with  white 
phosphorus . 

Conduct  additional  yaw  sonde  instrumented  flight  tests  of 
liquid-filled  projectiles. 

Perform  internal  flow  visualization  studies  in  laboratory. 

Investigate  experimentally,  the  effect  of  non-cylindrical 
containers  (i.e.,  cylinder.?  with  endcaps  that  are  ellipsoidal, 
conical,  etc.)  for  highly  viscous  liquids. 

Evaluate  longitudinal  baffles  to  reduce  destabilizing  moment 
due  to  highly  viscous  liquid-fills. 

Analyze  transient  effects  of  immiscible,  low  viscosity 
additives  in  reducing  instabilities  (yaw  sonde,  spin-up 
experiments,  theoretical,  etc.). 

d  Investigate  large  despin  moment  at  low  spin  rates  for 

viscoelastic  fluids. 
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Unsteady  and  gravity  effects  at  low  Reynolds  numbers. 

Non-linear,  unsteady  and  shape  effects  at  high  Reynolds 
numbers. 

Determine  Reynolds  number  limits  for  existing  codes. 

Establish  single  analytical  method  for  liquid-filled 
projectiles  that  handles  entire  Reynolds  number  range  (i. 
one  stop  method)  . 
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"WORKSHOP"  RECOMMENDATIONS 
(Not  Prioritized) 


1.  Obtain  more  detailed  effects  of  immiscible  additives  in 
preventing  flight  instabilities. 

2.  Study  effect  of  unconventional  geometry  (i.e.,  baffles, 
internal  payload  shape,  etc)  in  creating  and  preventing  flight 
instabilities . 

3.  Determine  influence  of  porous  media  in  reducing  flight 
instabilities . 

4.  Investigate  larger  range  of  visco-elastic  fluids  related  to 
flight  instabilities. 

5.  Understand  transient  effects  during  launch  and  flight. 

6.  Include  effect  of  chemical  reactions  during  flight. 

7.  Analyze  in-flight,  mixing  phenomena. 

8.  Evaluate  non-linear  dynamic  factors. 

9.  Evaluate  in-flight  thermal  effects. 

10.  obtain  additional  laboratory  and  flight  test  data  to 
provide  experimental  validation  of  theoretical  results. 

11.  Determine  influence  of  combined  internal  and  external  flov 
effects  on  flight  stability. 


Bl  ank 
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WORKSHOP  OH  PROBLEMS  OF  ROTATING  FLUIDS 

Program  Agenda 
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Monday,  April  22, 1991 

8:30  am  Welcome  Dr.  Donald  Austin 

Executive  Director,  AHPCRC 


8:45 


9:00 


9:30 


10:30 


Opening  Remarks 


Dr.  Daniel  Joseph 
Aerospace  Engineering  &  Mechanics,  U  of  MN 


Overview  of  CRDEC 
Research  Program  for  Liquid - 
Filled  Projectiles 

Analysis  and  Visualization  of  the  flow 
in  a  Spinning  and  Nutating  Container 
(Computer  Demonstration) 

Tour  of  Army  High  Performance  Computing 
Research  Center 


Mr.  Miles  Miller 
Chemical  Research  Development 
&  Engineering  Center 

Dr.  Thorwald  Herbert 
Ohio  State  University 


Dr.  George  Sell 
Director.  AHPCRC 


11:30  LUNCH 

1:00  pm  Tour  of  Fluid  Dynamics  Laboratory  Dr.  Daniel  Joseph 

Aerospace  Engineering  &  Mechanics.  U  of  MN 

2:30  A  Centrifugal  Spectrometer  Dr.  Harvey  Greenspan 

Massachusetts  Institute  of  Technology 

3:00  Moment  Exerted  by  a  Viscous  Liquid  In  a  Dr.  Charles  Murphy 

Spinning,  Coning  Container  Ballistics  Research  Laboratory 


Tuesday,  April  23. 1991 


8:30  am 

9:00 

9:30 

10:00 

10:30 

11:00 


Computational  Study  of  the  Unsteady  Flow  In  a  Dr.  Rlhua  LI 

Spinning  and  Nutating  Cylinder  Ohio  State  University 


Numerical  Simulations  of  Non-Cylindrlcal  Mr.  Michael  Nusca 

Liquid-Filled  Containers  Ballistics  Research  Laboratory 


Motion  of  Two  Immiscible  Fluids 
In  a  Spinning  and  Coning  Cylinder 

Direct  Measurement  of  Liquid 
Effects  Using  a  Moment  Balance 

Laboratory  Flight  Stability 
Evaluation  of  Production 
M825A1  Payload  Canisters 


Mr.  Mohamed  Selml 
Ohio  State  University 

Mr.  David  Hepner 
Ballistics  Research  Laboratory 

Mr.  John  Molnar 
Chemical  Research  Development 
&  Engineering  Center 


LUNCH 


(OVER) 
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1:00  pm 
1:30 

2:00 

2:30 

3:30 


Theory  and  Experiments  for 
Rotating  Porous  Media  Flow 

Effect  of  Interior  Canister  Wall  Roughness 
on  Liquid  Despin  Moment 


Instrumented  Flight  Tests 
Artillery  Projectiles  with 
Selected  Liquid-Fills 

Wrap-up  and  Final  Remarks 

Adjourn 


Dr.  Gene  Cooper 
Ballistics  Research  Laboratory 

Mr.  Daniel  Weber 
Chemical  Research  Development 
&  Engineering  Center 

Mr.  Miles  Miller 
Chemical  Research  Development 
&  Engineering  Center 

All 
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